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Goals 

Soil Moisture Data Assimilation 


Forecast Challenge 

- Available moisture affects humidity, 
sensible/latent heating, diurnal 
heating rate, and convection 

Objective 

• Improve soil moisture estimates 
for regional NWP applications and 
situational awareness 

- Improve LIS soil moisture by 
assimilating satellite retrievals 

- Use LIS output to initialize NWP 



Sensible HF Diff (LISMOD - CNTL) valid 080609/1600V013 


t'V': * *" 

Sensible HF 



2-m Dewp Diff (LISMOD - CNTL) valid 080609/1600V013 



Diff in PBL Hgt (LISMOD - CNTL) valid 080609/L600V013 CAPE Diff (LISMOD - CNTL) valid 080609/1600V013 


Impact of using high-resolution LIS boundary 
conditions in WRF (rather than NAM fields). 
From Case et at. 2008 
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Other Applications 
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SMOSand SMAP 

• L-band radiometers (and radars) can be used to estimate soil 

moisture in the top layer (~5 cm) of soil. 

-L-band (1.4 GHz) sees deeper in the soil, and performs better in dense 
vegetation than higher frequency instruments like AMSR-E (6-10 GHz). 

• Currently assimilating retrievals from European Space Agency's 
Soil Moisture and Ocean Salinity (SMOS) satellite 

• Preparing for assimilation of NASA Soil Moisture Active/Passive 
(SMAP) retrievals 

— Launch in early 2015 

-Combined (radar/radiometer) product available at a higher resolution (9 km) 

- We are members of the SMAP Early Adopters team 



Soil Moisture Instruments 
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Name 

AMSR-E 

SMOS 

SMAP 



Soil Moisture 
and Ocean 
Salinity 

Soil Moisture Active/Passive 



Agency 

NASA/JAXA 

ESA 

NASA 

Launch 


2009 

Nov. 2014 

Orbit 

Polar 

Polar 

Polar 

Sensor 

Type 

Passive 

Passive 

Passive 

Active 

Combined 

Frequency 

6.9 GHz 
(C-band) 

1.4 GHz 
(L-band) 

1.41GHz 

1.2 GHz 


Resolution 

56 km 

35-50 km 

36 km 

3 km 

9 km 

Accuracy 

6 cm 3 /cm 3 

4 cm 3 /cm 3 

4 cm 3 /cm 3 

6 cm 3 /cm 3 

4 cm 3 /cm 3 
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Soil Moisture Assimilation in LIS 


Noah LSM (in LIS) produces 
"operational" so I 
moisture/temperature states a 
surface fluxes of water and 
sensible/latent heat. 

Use data assimilation of 
SMOS/SMAP soil moisture 
retrievals to improve model 
states 

We plan to implement this in 
near-real-time SPoRT LIS to 
improve product for end users 


Land Information System 


r 

Topography, 



Soils 



(Static) 




J 


Vegetation, 

LAI 

(Monthly) 


Meteorological 
“Forcings’ 1 
(Hourly-3 hourly) 


Land Surface States 
(Variable) 


Land Surface Models (LSM) 
(15 min-1 hour) 

Noah, CLM, VIC 


w 



Surface Fluxes: 
Qnet,Qle,Qh,Qg 
Swnet, LWnet 


Soil States: 
Moisture, 
Temperature 
Profiles 


Water Balance: 
Runoff, Drainage, 
Baseflow (some) 


Surface States: 
Snow 

Skin Temperature 
LAI (some) 
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Data Assimilation with EnKF 

• LIS has a built in Ensemble Kalman Filter 

• EnKF combines the model background and observations to 
make analyses 

-Relative weighting is controlled by the specified observation error and by 
the ensemble spread 

• Implemented EnKF assimilation of SMOS data 

-Read SMOS data files 

-QC based on model state and data flags for precipitation, RFI, data quality, 
frozen soil, snow cover, and high vegetation 

-Ongoing experiments to tune run-time settings including perturbations, 
number of ensemble members 

-Bias correction by CDF Matching 

-Capability of implementing landcover-dependent correction. 



8 


Bias Correction 


• Initial tests had large dry bias in observations, so that only 
extreme rain events had correct sign. 



Discussions with other researchers confirmed need for bias 


correction 

Innovations 


Increments 


Implemented CDF matching 
correction for SMOS 
retrievals. 

Assimilating retrievals (not 
radiances) lets us use 
established methodology 


Uncorrected innovations (observations minus model) and 
increments. Red=dry bias in retrievals. 



9 


Bias Correction 


CDFs of Soil Moisture Observations and Model Correction Curves for 3 vegetation categories 


CDFs 



Modeled vs. Observed Soil Moisture 



* LIS can apply point-by-point correction curves. To increase the background dataset 
size, we are aggregating points by landcover type. We will also explore correction 
at each point and aggregating by soil type. 
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Irrigation Case Study 


0.09 0.12 0.15 0.18 0.21 0.2+ 0.27 



Model soil moisture 
concentration forced only 
by precipitation and 
misses magnitude of 
irrigation-saturated MS 
Valley 


SMOS observes 
irrigated fields 


Blended analysis of model 
and observations better 
represent irrigated area 
and should result in 
improved weather and 
hydrologic modeling 



From U. of Frankfurt-FAO (Ozdogan 2001) 


Test Impact on NWP 
using coupled LIS-WRF 

Implications for regional 
cl mate model ng 

— Impacts of changing land- 
use, precipitation patterns 
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Current and future plans 

• Testing DA for higher resolutions (grid cells « observation size) 

-Use in SPoRT high-resolution LIS runs. 

• Validate analyses 

-TAMU North American Soil Moisture Database 

• Test Impact of assimilating SMOS retrievals on NWP using 
coupled runs in NU-WRF 

-Impact on boundary layer for a quiescent day 
-Active convection case 
-Validation over a longer time period 
-Look at both sensitivity and forecast accuracy 

• Implement with SMAP retrievals 

Questions? 
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Perturbation Tests 



Background 12Z 1 Apr 2013 (pert= .01) 


0.15 0.18 0.21 



Innovotion 12Z 1 Apr 2013 (pert=.01) 


Increment (pert=.01) Anolysis 12Z 1 Apr 2013 (pert=.01) 




Bockground 12Z 1 Apr 2013 (pert=.02) 



Innovotion 12Z 1 Apr 2013 (pert=.02) 


-0.15 -0,1 -0.05 -0.01 0.01 



Increment (pert= .02) 





Analysis 12Z 1 Apr 2013 (pert=.02) 


Background 12Z 1 Apr 2013 (pert=.04) Innovotion 12Z 1 Apr 2013 (pert= .04) Increment (pert=.04) Analysis 12Z 1 Apr 2013 (pert=.04) 





